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The concept of effective particles as degrees of freedom in a relativistic 
quantum field theory is defined using a non-perturbative renormalization 
group procedure for Hamiltonians. However, every candidate for a basic 
physical theory appears to require an initial perturbative search for the set 
of interaction terms that may provide a basis with which the full effective 
theory Hamiltonian could be constructed in a series of successive approx- 
imations. This article describes the required perturbative expansion and 
illustrates it with a set of general 4th-order formulae. 



1. Introduction 

The construction of relativistic interactions of effective particles that is 
described in this article is designed for a quantum theory in which one can 
introduce a front form (FF) of Hamiltonian dynamics [T] . The construction 
draws on the previous work on similarity renormalization group (SRG) pro- 
cedure for Hamiltonians [2], flow equation [3], and renormalization group 
procedure for effective particles (RGPEP) [1]. It is the latter that is used 
here. 

The key limitation of the RGPEP in comparison to the general SRG 
procedure p] is that in RGPEP one focuses on coefficients of operators in 
an operator basis built from products of creation and annihilation operators 
for effective particle^ rather than on the Hamiltonian matrix elements. The 

^ For brevity, the creation and annihilation operators will be sometimes commonly 
called particle operators. 



(1) 
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requirement that the effective dynamics is obtained by universally rotating 
the bare particle operators to corresponding effective ones imposes unitary 
constraints on the effective theory. These constraints are not necessarily 
satisfied when one defines an effective Hamiltonian by rotating its matrix 
in some basis. In the matrix, many operators may contribute to one and 
the same matrix element and a rotation of a matrix does not have to be 
equivalent to a rotation of particle operators. In the effective particle the- 
ory, the matrix rotations one considers result only from rotating particle 
operators. This simplification allows one to encode many matrix elements 
in a small number of operators. In practice, the simplification begins to 
manifest itself when one considers states with more particles than 2 or 3. 
For example, a two-body effective potential term involves a product of only 
four fields and a function, say v, of their four arguments, such as in Sees. II 
E and X C in [5j. In contrast, a matrix element of the same term between 
states of many particles involves many copies of the function v with different 
arguments. In addition, since the field operators involve both creation and 
annihilation operators, even an operator built from just four field operators 
requires a large space of states to study its matrix elements and uncover its 
actual operator structure result. 

Derivation of the perturbative formulae for relativistic interactions of ef- 
fective particles starts from the non-perturbative RGPEP equation that has 
been mentioned before in [S] on the basis of jTJ. The starting RGPEP equa- 
tion resembles Wegner's j3] (for reviews, see [51[^) in its double commutator 
structure. Somewhat different perturbative expansions could be obtained 
in the RGPEP starting from the non-perturbative equations that involve 
multiple commutators. Such equations are described in Appendix [A} 

The relativistic nature of effective interactions is achieved using RGPEP 
by respecting the 7 kinematical symmetries of the FF of Hamiltonian dy- 
namics [Ij and cluster properties |10j in quantum field theory. The starting 
non-perturbative equation of RGPEP is designed to preserve these features. 
They are also preserved in the perturbative expansion. The running cutoff 
parameter of RGPEP limits only changes of the invariant mass of inter- 
acting particles. These features are required, for example, in application 
to QCD where one desires to simultaneously explain the constituent quark 
model classification of hadrons pjj and their quite different picture in the 
parton model |12j . 

Section [2] provides a brief account of the RGPEP and Section [3] describes 
the perturbative solution up to the 4th order. Section [4] concludes the paper. 
Appendix |A] describes multi-commutator fiow equations and Appendix [B] 
provides universal formulae for coefficients in the perturbative expansion of 
effective Hamiltonian interaction terms. 
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2. Non-perturbative RGPEP [7J 

Consider a local field theory in which fields on a given hyper-plane in 
the Minkowski space-time are expanded into their Fourier components. The 
components have interpretation of creation or annihilation operators of the 
field quanta. These quanta are called the bare, or initial particles, and 
the operators that create or annihilate them are called bare particle oper- 
ators. The bare particle operators are generically denoted by go- A one- 
bare-particle state is obtained by acting with a creation operator go on the 
vacuum state |0). The bare particle operators are used to build a basis in 
the Fock-space by acting with products of creation operators go on |0). 

Hamiltonian densities of canonical quantum field theories are built from 
products of fields and their derivatives. A Hamiltonian obtained by inte- 
grating such a density over a space-time hyperplane is a combination of 
products of operators go with coefficients cq. If a product contains n par- 
ticle operators, the coefficient has n arguments. Each of the arguments 
represents a complete set of quantum numbers carried by a corresponding 
particle. One sums or integrates over these arguments in the Hamiltonian 
terms. 

RGPEP introduces effective particle operators 

g. = UsQ^uU (1) 

labeled by the parameter s which plays the role of a renormalization group 
parameter. This means that s labels a family of Hamiltonians that all cor- 
respond to one and the same theory but are expressed in terms of differently 
defined degrees of freedom. All kinematical quantum numbers of g are the 
same on both sides of Eq. ([l]), i.e., irrespective of the value of s. 

The parameter s has dimension of length and ranges from to any finite 
number of choice, including arbitrarily large numbers. The larger s the 
harder calculation of the corresponding Hamiltonian. s — t- oo corresponds 
to solving for the spectrum of a theory. This is exemplified in the elementary 
model discussed below. In complex theories, where there is little hope for 
obtaining analytic solutions and one has to use numerical methods, the 
parameter s may hopefully be kept, on the one hand, sufficiently small 
to maintain an analytically controlled connection with the quantum field 
theory one starts from, and, on the other hand, sufficiently large to enable 
numerical calculations using computers. 

Physically, s has the interpretation of the characteristic size of the ef- 
fective particles. The interpretation follows from the fact that effective 
interactions contain form factors that limit how far off energy shell the in- 
teractions can extendj^ The width of the form factors is determined by 1/s. 

The energy scale is introduced using the spectrum of a free part of the Hamiltonian. 
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Consequently, s = corresponds to the concept of point-like, bare particles 
in the case of a local quantum field theor}|^ The effective Hamiltonian that 
corresponds to scale s is band-diagonal on the energy scale (more precisely, 
on the scale of invariant masses of particles involved in an interaction) and 
the band width is ~ The principle of using the band-diagonal structure 
for the purpose of renormalization is formulated in [2]. 

Let the initial Hamiltonian including counterterm^ be 7io{qo)- For di- 
mensional and notational reasons, it is convenient to use the parameter 
t = and label Hamiltonians and other operators with t rather than s 
itself. The RGPEP demands that 

nt{qt) = -Hoiqo) . (2) 

This means that one changes particle operators from go to qt and at the 
same time the coefficients cq are changed to ct so that the Hamiltonian 
operator stays unchanged. Consequently, the wave functions of its eigen- 
states, in a Fock-space basis built using qt, depend on t. This design serves 
the purpose of simultaneously explaining the canonical picture of Hamilto- 
nian eigenstates as built from bare, point-like particles, and the constituent 
picture in terms of extended, effective particles. Such setup is desired for 
solving QCD and comparing solutions with experimental data. 

2.1. Generator 
By differentiating both sides of 

Utiao) = ulno{ao)l(t, (3) 

with respect to t, which is denoted by a prime, one obtains 

n'tiao) = [-UlulUtiao)]. (4) 

The product 

Qt = -uH (5) 

is called a generator. The RGPEP generator is defined by the formula 

Gt = [nf,npt], (6) 

When the RGPEP is apphed to a non-local effective theory, the starting value of 
s = corresponds to the initial size of non-locality. 

The counterterms have to be derived and the derivation is a part of the RGPEP but 
the initial Hamiltonian is meant to contain the required counterterms, i.e., the ones 
to be found using the RGPEP. 
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whose right-hand side will be explained below. The definition allows for 
adding arbitrary multiples of 7it{ao) or other operators that commute with 
l-Lt{o-o) but such additions to Qt are ignored as immaterial here. 

Since Qt is defined as a commutator of two formally hermitian operators, 
it is formally anti-hermitian, and lit is formally unitary. The word "formal" 
is used here because precise definitions require regularization (see the quoted 
literature). 

In what follows below, if an operator is expressed in terms of the basis 
operators qq, they are no longer indicated as its arguments. However, one 
should remember that the effective Hamiltonian is obtained by making the 
substitution, H = litiqt) = T~it{<lQ — ^ In the notation for operators go 
themselves, the subscript is also omitted, so that g = go- 



2.2. Uf in Eq. ^ 

The Hamiltonian T-Lf^ called the free Hamiltonian, is the part of a full 
Hamiltonian at t = = that has the form 

'Hf = Y.P^qW (7) 



and does not depend on the interactions. More precisely, T~Lf is defined as 
the part of % that one is left with when the coupling constants in the initial 
theory are set to and only the terms that are bilinear in fields are kept. 
The sum extends over all particle species and their quantum numbers so 
that the sum also includes integration over particle momenta p. The minus 
component of a particle momentum, = Pi — pf, in Eq- P^, 



P7 = 2t!±^ . (8) 

pT 

is the FF energy that corresponds to a free particle with the mass nii and 
kinematical momentum components pf and pj~ . Typically, different species 
of particles have different masses. 



2.3. Upt in Eq. 

The operator T-Lpt is uniquely defined once the Hamiltonian 1-Lt is spec- 
ified. If the Hamiltonian T-Lt is of the form 

oo 

n=2 ii,i2,...,i„ 
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where the coefficients ct{ii, ■■■,in) are to be found using RGPEP, the oper- 
ator Tipt is defined by 

OO / n \ ^ 

= E E '^t{^l,.■■,in) [lY^Pt] li-'-Q^n- (10) 

n=2 11,12,..., in \ k=l / 

Thus, Upt differs from Tit by multipHcation of its terms by the square 
of +- component of total momentum carried by the particles in a term. 
This kinematical momentum is specified by the operator content of a term 
irrespective of the value of RGPEP parameter t = s^. 

The definition of Upt secures that both sides of Eq. Q, which now 
reads 

j^-Ht = mf,npt],'Ht] , (11) 

behave in the same way with respect to operations of kinematical LF sym- 
metries. This implies that the effective particle size parameter s is invariant 
with respect to the FF kinematical subgroup of the Poincare group. 

2.4- Construction of counterterms 



The RGPEP Eq. (11) predicts the coefficients q of products of qt in 
Utiqt) provided the initial condition for 7io{Qo) is available. However, when 
the initial Hamiltonian involves divergences, such as the ones due to a dy- 
namical coupling of infinitely many degrees of freedom over an infinite range 
of scales, it cannot be considered a good initial condition. One can cut the 
infinities off by regularization, limiting the space of states and/or limiting 
the range of scales involved in interaction. With the cutoffs, the solutions 
lead to effective Hamiltonians that involve huge numbers (or zeros) instead 
of meaningful terms. The huge numbers result from the ratios of cutoffs to 
physical parameters. Zeros result from inverses of the huge ratios. 

The situation is understood here as the one in which the initial Hamil- 
tonian may reflect a physically interesting structure that is relevant when 
one imposes so small cutoffs on the dynamics that the divergences are re- 
placed by relatively small terms and these terms do not obscure the initial 
structure beyond recognition. The question is then how a good effective 
Hamiltonian should depend on the cutoffs so that its predictions do not. 
If the predictions did depend on the cutoffs, one could not treat the cutoff 
theory as equivalent to some fundamental one. By the word "fundamental" 
it is meant here that the range of validity of a "fundamental" theory is 
expected to be vastly greater than the range limited by adjustable cutoffs. 

Following \13\ [TS] and [2] , the RGPEP involves determination of addi- 
tional terms, called counterterms, that need to be included in the initial 
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Hamiltonian in order to render the effective theory which may correspond 
to the Hamiltonian initiaUy suggested as vahd. The construction of coun- 
terterms follows the rules of SRG procedure [21 01 [7] . There is no explicit 
condition of widening of the band toward high energies (large invariant 
masses), but the narrowness of the Hamiltonian (see next section) enables 
one to search for counterterms on the basis of a condition that all matrix 
elements of a Hamiltonian corresponding to some finite size s of effective 
particles between states of finite invariant masses, are independent of the ini- 
tial cutoffs. This condition provides as many equations to solve as there are 
matrix elements to consider. The number of equations may be large enough 
to specify the structure of counterterms up to a limited set of unknown fi- 
nite numbers or functions. These numbers or functions, called finite parts of 
counterterms, can be constrained by symmetries [15] and phenomenology. 



In general, the only method for solving Eq. (11) is to use successive 

(n) 

approximations. One hopes that a candidate number n for a solution, T-Ll , 
can be inserted on the right-hand side to render a new candidate ^j""*^^^ on 
the left-hand side, and ^^"^^^ can be subsequently inserted on the right- 
hand side in place of Til"^ to render T^^"^^^ on the left-hand side, and so on. 
More precisely, the hope is that a well-defined solution is approximated with 
increasing accuracy when n increases. Of course, one has to study theories 
case-by-case in order to establish if the RGPEP sequence converges [16j. In 
principle, one may obtain not only solutions for Tit with limit cycles |17j . 
instead of just fixed points, but also the solutions with chaotic behavior as 
functions of t ^HHj. 

The conceptual and computational difficulty is that in each step of con- 
structing successive approximations one may have to modify the countert- 
erms and the modification is not dictated by the iteration procedure alone. 
Namely, one has to inspect dependence of the small-invariant-mass matrix 
elements of H^"^ at large t on regularization present in T^j"^ at t = and 

(n) 

attempt to counter this dependence by modifying the Til at t = 0. The 
modifications may amount to redefinitions of a finite number of coupling 
constants or functions of particle quantum numbers. However, when one 
confronts the issue of relativistic description of confinement, the question of 
convergence is not answered in any form yet and the mechanism by which 
effective Hamiltonians develop confining forces at finite t is strictly speaking 
not known. 

The attractive feature of the RGPEP equation is that it allows for split- 
ting of the integration process into many small steps, accompanied with 
rescaling of invariant masses so that one always operates with dimension- 
less quantities sAi, where is a free invariant mass of the interacting 
effective particles and s is their size parameter. Each small step can be 
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executed using variety of well-known mathematical methods as in the case 
of the original Wilsonian procedure \13\ I14j . Thus, the RGPEP provides 
a new tool for studying universality in relativistic quantum field theories 
(see Appendix A in [7J). This feature is of interest because a well-defined 
effective QCD Hamiltonian may be sought using the RGPEP irrespective 
of many details in the initial Hamiltonian. Namely, it should be sufficient 
to start with any Hamiltonian in a suitable universality class. 

Nevertheless, in order to begin the process of constructing solutions of 
the RGPEP equations as outlined above, one has to suggest the terms that 
should be included in Ti^^^ and can be expected to have significant coeffi- 

(n) 

cients in T-L^ with large n. Initially, the only available tool for gathering 
such information about important terms is perturbation theory. This arti- 
cle illustrates a general algorithm for generating perturbative formulae for 
Ht- The illustration includes a set of perturbative formula up to 4th order, 
derived in Section [3l 



2. 5. Narrowness 



Solving Eq. (11) by successive approximations involves restrictions on 
how many terms are kept in the Hamiltonian. Each of these terms is written 
as a product of quantum fields. A priori, a FF Tit contains infinitely many 
terms [5j. A good approximate operator solution must be little sensitive to 
the terms that are missing. Assuming that one can identify a set of dominant 



operators, the mechanism of narrowing the Hamiltonians through Eq. (11) 



can be seen using an equation that results from projecting Eq. (11) on a 
subspace of the Fock space. 

Following Appendix C in Ref. [7j, one can introduce a projector on 
a subspace of physical interest in the Fock space. Both the subspace and 
projector are denoted by R. The projection involves fixing total kinematical 
momentum of states and a cutoff on their total invariant mass so that the 
concept of a trace becomes well-defined. Let 

n = RUtR. (12) 

The correspondingly projected RGPEP equation reads 

n' = [[Uf^npin] . (13) 

The condition that the trace of Ti^ does not depend on t, implies for matrix 
elements T-lmn = {m\T-l\n) in the basis built from eigenstates \m) oiHf with 
eigenvalues PJ^-, that 



. mn 




Y^mrrr^n?) = -Y.^Pju - Pjm?\^Ikm?2P+l < , (14) 
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where P^^^ is the total +-momentum of the particles that are involved in 
the interaction. In terms of the invariant masses, 

= -2Y,{Ml^-Ml,^f\-Hikm?<^. (15) 

km 

where Mkm denotes an invariant mass of the particles in state labeled with 
k that are transformed through the interaction T-Li to particles in the state 
labeled by m. Change of order of subscripts results in an invariant mass 
of particles in a different state. Spectators do not contribute to the invari- 
ant masses. The calculations are explicitly invariant under 7 kinematical 
transformations of the FF of Hamiltonian dynamics. 



Eq. (15) means that the sum of moduli squared of all matrix elements 
of the interaction Hamiltonian decreases as t increases until all off-diagonal 
matrix elements of the interaction Hamiltonian between states with different 
free invariant masses vanish. Section [3] shows this feature in a perturbative 
expansion for Ht in powers of Hit- 

2. 6. Transformation lAt 
Transformation Ut is a solution to 

u[ = -Ut [nf,npt\. (16) 

The initial condition isUq = 1 . Successive approximations generate Ut from 
a solution for Tit- The general solution is 

Ut = Texp(^- J^dT[nf,npr]y (17) 

where T orders operators from left to right in the order from a smallest to a 



largest of their arguments t. Perturbative expansion of Eq. (17) provides a 
perturbative solution for Ut- Knowing Ut, one can calculate operators that 
create or annihilate effective particles of size s using Eq. M with t = s^. 



3. Perturbative formulae 



Solution of Eq. (11) for coefficients of products of effective particle 
operators in Tit, is described in this Section using the notation that proved 
useful befor^ in which Eq. ( [ll| takes the form 



'^'t ab = -ob^ 'Hit ab + X] ^P"'"' "'^ + P'"' ^-f* ^-f* ■ 



(18) 



E.g., see (20J, Sec. II.B, or [2T], Sec. III.C. 
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The letters a, b, and x, denote configurations of particles. A configuration 
is a collection of quantum numbers that label particle operators. 

An example illustrating the concept of configuration is shown in Fig. [TJ 
In Fig. [T| the configuration a includes particle operator labels contained in 
the sets denoted by 1 and 2, the configuration x includes particle operator 
labels contained in the sets denoted by 1, 5 and 4, and the configuration 
b includes particle operator labels contained in the sets denoted by 3 and 
4. A subscript such as ax refers to a coefficient of the product of particle 
operators in a Hamiltonian term Tij that changes a configuration x to a 
configuration a, etc. In Fig. [T} Hit ax would correspond to a term in Hit 
that involves a product of creation operators labeled with quantum numbers 
in set 2 and a product of annihilation operators with labels contained in 
sets 5 and 4. Tijtxb would correspond to a term that involves a product 
of creation operators labeled with quantum numbers in sets 1 and 5 and a 
product of annihilation operators with labels contained in set 3. Symbol ax 
as a term in the equation denotes a difference of invariant masses squared, 
ax = M-ax ~ -^xa- -M-ax denotes an invariant mass corresponding to only 
those particle operator labels in a configuration a that emerge as a result 
of interaction from the configuration x. Thus, ax = —xa. For example, 
in Fig. [T| A^aa; denotes the invariant mass corresponding to the particle 
operator labels in set 2 and A4xa denotes the invariant mass corresponding 
to the particle operator labels in sets 5 and 4. Spectators do not contribute 
to these invariant masses. Thus, in Fig. [T| Aiax and Aixa do not include 
particle operator labels from set 1, while Aixb and Aitx do not include 
particle operator labels from set 4. The invariant masses are evaluated 
using kinematical momenta and eigenvalues oi Tif. A symbol such as pax 
denotes the sum of + components of all momenta that appear in the labels 
of creation operators (or of all annihilation operators, which is the same due 
to the momentum conservation) in a product in the Hamiltonian term that 
changes the configuration x to a, etc. In Fig. [T| pax denotes the sum of + 
components of momenta contained in set 2, which is the same as the sum 
of + components of momenta contained in sets 5 and 4 together, while p^t 
denotes the sum of + components of momenta contained in sets 1 and 5, 
which is the same as the sum of + components of momenta contained in 
set 3. 



Key features of Eq. (18) are following. Since Eq. (18) only refers to 
quantum numbers that label particle operators, the Hamiltonian operators 
that solve Eq. (18) act in the entire Fock space. This is how the RGPEP 
procedure avoids the limitation to an a priori limited set of states that one 
is forced to work with using equations for Hamiltonian matrices. There 
are no disconnected terms in T-Lt, since the second term in Eq. (18) results 
from a commutator. The commutator also implies that one keeps only the 
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Fig. 1. Example of the RGPEP operator calculus, see Eq. ( 18 1 and its description. 



terms that result from commuting at least one annihilation operator with 
one creation operator in the product Hit ax T~{-itxb- Also, Tit does not explic- 
itly depend on the eigenvalues of a full Hamiltonian. Such dependence is 
a serious limitation of all procedures that are in principle based on Gaus- 
sian elimination, or, in the path integral formulation, on integrating out 
the states outside a cut off in the space of states [T^. Namely, the larger 
an eigenvalue the worse accuracy of the procedure for integrating out the 
components above a cutoff. 

The problem of integrating out states above a cutoff is circumvented 
in the RGPEP at the level of constructing counterterms and evaluating 
effective Hamiltonians. A finite remnant of this problem remains in an 
eigenvalue equation for Titiqt), when one limits the space of states of effec- 
tive particles that are included in a non-perturbative diagonalization using 
computers and when one fits finite parts of counterterms by comparison 
of theoretical results with experiment. This residual issue must be dealt 
with in the diagonalization procedure and typically should not cause prob- 
lems because interactions do not grow with energy fast enough to generate 
any divergences. If they did cause divergences, a separate search for fixed 
points or other in principle possible behavior would have to be arranged 
(see Appendix B in 0). 



3. 1 . Form factors 



Eq. (18) suggests that one writes 

T'I'tab = 

where 

ft ab 



—tab"^ (-1 

e ytab-, 



-tab"^ 



(19) 



(20) 



is called a form factor. The width of the form factor as a function of the 
invariant mass change is 1/s. Therefore, the parameter s has the inter- 
pretation of a size of effective particles. It determines how far off shell an 
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interaction can reach in terms of the invariant mass. The form factor se- 
cures narrowness of the effective Hamiltonian order-by-order in perturbation 
theory. 



Eq. (18) imphes that Qtab satisfies the equation 

^iab = ''^^taxbGtaxGtxb , (21) 

X 

where 

Ataxb = {Pax ax+p,, bx) ^-t{ax^+-b^-ab^) _ (22) 



Eq. (21) is solved below order-by-order in a perturbative expansion using 



powers of a single coupling constant g. 

3. 2. Expansion 

When one expands 

Qt = nf + ggti + g^gt2 + g^Gt3 + ..-, (23) 

the perturbative formulae for calculating Gtu is obtained in the form 

n-l 

S'tnab = ''^^''^^^taxbStkaxGtn-kxb- (24) 
fc=l X 

It yields order-by-order the following differential equations 

Q'tlab = 0, (25) 
Gt2ab = ''^^taxbGtlaxGtlxb, (26) 

X 

G'tSab = At axb {Gtl ax Gt2 xb + Gt2 ax Gtl xb) , (27) 

X 

G'u ab — A-t axb {Gtl ax Gt3 xb + Gt2 ax Gt2 xb + GtZ axGtlxb) ■ (28) 

X 

The four terms indicate a generic pattern. This pattern is different from the 
patterns described in perturbative expansions studied using matrix models 
in Refs. [23] and [23]. The main reasons for the patterns to differ are the 
constancy ofHf, which is of order 1, and the definition oiTip in Eq. (|11[ ), 
which leads to the factors of -|— momentum of active particles in Eq. (\22^ . 

The same factors of p'^ appear in the procedure studied in Ref. j24|. 
but there they are accompanied by appearance of the derivative of G' on 



SGlazekRGPEP20120420 printed on January 30, 2013 



13 



the right-hand side of the RGPEP equation for %' ^ in the generator. The 
derivatives in the generator cause additional difficulties in solving the RG- 
PEP equation non-perturbatively. Here, the origin of factors of is in 
T-Lpt- The absence of derivatives in the generator remedies the additional 
difficulties encountered in Ref. [24J. 



General formulae for solutions of Eq. (24) are listed below order-by- 
order for all terms up to the 4th order, i.e., up to 2 or 3 loops in theories of 
the type 0^ or (p^, respectively, where the power of ^ indicates how many 
field operators can be present in a term. This covers in principle all cases of 
interest in physics (except for gravity, because of the FF limitation to the 
Minkowski metric). For example, the expansion to 4th order is in principle 
sufficient to study interactions of constituent quarks in hadrons including 
the leading effects of self-interaction, gluon exchange, and running coupling 
constant. Such studies should clarify if the 4th order RGPEP is capable of 
generating information about the essentially non-canonical effective QCD 
interaction terms discussed in Ref. or one must explore orders higher 
than 4 to see the new terms. Terms of higher orders can be generated as 
required according to the pattern visible in what follows. 

3.3. Solutions, order-by- order 



Appendix B describes details of solving Eqs. (25) to (28). This section lists 



the results, introducing the elements of notation developed in Appendix [B| 
3.3.1. Ist-order solution 

The first-order solution does not depend on t, 



Hlab 



Goi ab- (29) 



The subscript refers to t = 0. Operators with subscript are contained 
in the initial Hamiltonian. They provide the initial conditions. 

The first-order initial condition includes a bare coupling constant, as a 
coefficient in front of a definite canonical operator. Counterterms of higher 
orders typically include the same operator. The net result of including 
counterterms is a change of value of g. The magnitude of change depends 
on the interaction terms included in and regularization method adopted for 
the initial canonical Hamiltonian. 

3.3.2. 2nd-order solution 

The second-order solution reads 

Gt2ab = G 02 ab + '^2 ^twrb^^ ^02 axb ■ (30) 
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The initial condition Qo2ab typically includes instantaneous interactions, 
such as a FF counterpart of the Coulomb interaction term in the standard 
form of dynamics, or instantaneous fermion terms that are unique to the FF 
of dynamics, due to constraints. The RGPEP may provide the higher-order 
terms that have similar structure but whose coefficients are not necessarily 
constrained according to the simple Heisenberg equations of motion that 
one might expect from the analogy with classical field equations [22]. For 
example, quantum interactions may lead to anomalies. The initial condition 
Go2 ab also includes self-interaction counterterms. 

The operator structure Gmaxb is defined in Eq. ( B.6 ). It results from 
action of first-order terms twice. The coefficient B'^^^^ is defined in Eq. 
(B.12). The superscript convention is explained in Appendix B.l It is well- 



known that the second-order terms reproduce standard second-order results 
for observables to the extent they are available in all theories relevant to 
physics. 



3.3.3. 3rd-order solution 

The third-order solution reads 



H2, ab 



+ 



E' 

xy 



(123,0) ^ 
t axb ^03 axb 



„{124,(234,0)) „(134,{123,0)) 
t axyb ' t axyb 



^03 axyb • 



(31) 



The term ^03a6 is an initial condition for an operator that does not appear 
in a canonical Hamiltonian. Namely, it is a third-order coupling constant 
counterterm. Operators Q^zaxb involve a product of two operators from 
the initial Hamiltonian: one of first order and another one of second order, 
see Eq. (B.7). The last operator, Gozaxyb-, is a product of three first-order 



operators, see Eq. 



R8|). The coefficients sj^.f and ^r^^T^""^ are 



,(134,(123,0)) 



defined in Eqs. ( B.lSp and (B.14) in Appendix B.l 



The third-order solution of Eq. (31 ) slightly differs from the third-order 



formula previously studied using RGPEP in the cases of scalar theories in 6 
dimensions [20] and in QCD in 4 dimensions [21]. Although the previous ex- 
perience suggests that Eq. (31) leads to the same result concerning asymp- 



totic freedom, the finite invariant-mass details of Eq. (31) are of special 



interest since they contribute to the fourth-order terms. The fourth-order 
terms required prohibitively complex equations in the previous applications 
of RGPEP. A relatively compact fourth-order solution is now made available 
in the next section. 
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3.3.4. 4th-order solution 

The fourth-order solution reads 

,(123,0) 



Hiab 



+ 

xy 

+ 

xyz 

+ B, 



X 

(124,(234,0)) „(134,(123,0)) 



t axyb 



GoAab + ^ Bl'^xb"' QoAaxb 

\ ^ r R(j-^'±,l,^o'i,'J;; , p 
/ . axyb ^ 
xy 



^04 



axyb 



„(125, (235, (345,0))) r^vJ-^^A^'^'JA^J'*!'-';;; i tdV 
taxyzb ' taxyzb ' taxyzb 

(145,(134,(123,0))) n(135, (123,0), (345,0)) 



t axyzb 



t axyzb 



(145,(124,(234,0))) 
^04 axyzb • 

(32) 



The initial condition ^04ab is the fourth-order counterterm. It includes the 
self-interaction and coupling- const ant counterterms that involve additional 
particles in an interaction. It also includes box-diagram counterterms, where 
they are needed. The operator Q^Aaxb is defined in Eq. (B.9), the operator 



QoA axyb is defined in Eq. (B.IO), and the operator GoAaxyzb is defined in 
Eq.(B.ll). The new coefficients in fourth order terms are defined in Eqs. 



(B.15) to (B.19) 



The general fourth-order result opens a door to many specific studies. 
Perhaps one of the most instructive ones would be a calculation of T family 
masses and decay width, which is likely to shed some light on the dynamics 
of gluons in heavy quarkonia, cf. f25|. 



4. Conclusion 

The perturbative RGPEP formulae of Section [3] allow one to study low- 
order divergences and effective interactions in specific theories with infinites- 
imal coupling constants and extreme cutoffs in their canonical Hamiltonians. 
They also provide a tool for testing the formal RGPEP claim of narrowness 
of effective theories described in Section [2j Such tests are essential because 
of anomalies, i.e., the new Hamiltonian terms that do not correspond to 
classical Lagrangians and instead result from imposing some regularization 
on a formal quantum theory. One potential source of anomalies that needs 
to be studied is a lower bound on +-momenta of particles in a Fourier anal- 
ysis of a quantum field. Strictly speaking, such lower bound violates boost 
invariance of a formal theory. 

Regarding systematic checks of existence of anomalies in the FF of 
Hamiltonian dynamics, there is currently no alternative known to the au- 
thor to inspecting theories order-by-order using the RGPEP and finding out 
term-by-term if and how soon a specific regularization method in a given 
theory may generate anomalies in the weak-coupling expansion. Moreover, 
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there is currently nothing known yet for certain about the anomahes gener- 
ated in the FF of Hamiltonian dynamics using RGPEP. 

The subject certainly deserves a study in view of the fact that the FF 
vacuum problem is differently formulated from the instant form vacuum 
problem. For example, if one imposes the above mentioned lower bound on 
+- momenta of particles in a Fourier analysis of quantum fields, the mathe- 
matical vacuum state, i.e., the state that is annihilated by bare annihilation 
operators, is an exact eigenstate of the full Hamiltonian with eigenvalue 0. 
If this state can play the role of the true ground state, the mechanisms of 
symmetry breaking and mass generation in effective theories may in prin- 
ciple turn out not to be associated with the vacuum expectation values of 
normal-ordered products of quantum fields but with some new terms in the 
Hamiltonians [5]. The perturbative RGPEP formulae for relativistic inter- 



actions of effective particles derived in Section 3.3 provide a new tool for 
studies of such hypotheses. 

Most interestingly, however, the perturbative studies may help in iden- 
tifying a finite set of operators that one can use to approximately solve the 



RGPEP Eq. ( 18 ) non-perturbatively, i.e., in terms of the coefficients in front 
of the identified operators. As functions of t, these coefficients would be ex- 
pected to evolve from a set of constants and regulating functions at t = 
to a set of finite, non-trivial functions of particle quantum-numbers at some 
to = Sq, where sq denotes some size of the effective constituents that are 
most suitable as degrees of freedom for specific purposes of phenomenology. 

Appendix A 

Multicommutator RGPEP 

One can obtain narrow relativistic Hamiltonians using equations with 
an odd number of commutators. Namely, one can write 

j^nt = [[nf,[.:nf,[nf,npt]...]],nt] , (a.i) 

including 2n — 1 operators Tif in a sequence of commutators with n > 1 
and t = s^". Correspondingly, one can define the operator Tipt by writing 

OO / Tl \ 

^* = E E ct{^u■.■,^n) [lYpt] i---'?^.. (A-2) 

n=2 ii,i2,...,in \ fc=l / 

instead of Eq. (10). The resulting counterpart of Eq. (18) reads (cf. |7j, 
footnote 43) 
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The narrowness is obtained for the matrices of projected Hamiltonians as 



in Section 2.5 However, instead of condition ( |15[) , one obtains a similar 
inequahty with {Ml^ - TW^^)^" in place of - Ml^^f . The per- 

turbative expansion proceeds without any qualitative change, but the form 
factors in Section [3.11 are of the form 



fi 



t ab 



(A.4) 



Appendix B 

Calculation of the ^th- order terms 

This Appendix contains a derivation of the fourth-order term described 



in Section 3.3 It also explains the pertinent notation. The calculation 



involves all terms of orders lower than 4. Using solutions given in Eqs. (29) 
(30), and (pll), one obtains from Eq. (28) that 



Uab 



t axb {yOlax y03xb 



+ ^ axb Bf^^yf^^^ {Golax Go3 xyb + Go2 ax Go2 xyb) 



+ ^03 ax ^01x6 + Gmax Go2xb) 



xy 



+ ^ At ayb -Bf axy (^03 axy GQlyb + Go2 axy Go2 yb) 
xy 

^ A [0(124,(234,0)) „(134,(123,0))1 ^ ^ 

-t- 2_^^taxb Y^txyzb ^ ^txyzb \ f 01 ax y03 xj/zb 

xyz 

^ Sr A [0(124,(234,0)) „(134,(123,0))1 ^ ^ 

+ 2^^tazb yi^taxyz + ^taxyz \ ^03 axj/z t/Ol zb 

xyz 

^ A r(123,0) „(123,0) ^ ^ 

+ Z^^taybl^taxy ^tyzb axy y 02 yzb , 



xyz 



where 



Go2axb — GoiaxGoixb , 
y03 axb ^01 ax ^02 xb + ^02 ax Goi xb ) 
^03 axyb = Goi ax Goi xy Goi yb ■ 



(B.l) 



(B.2) 
(B.3) 
(B.4) 



Inserting these definitions and grouping coefficients in front of the same 
operators, one obtains the derivative of fourth-order terms in the form 

^t4a6 = At {Goi ax ^03 xb + ^03 ax Gqi xb + ^02 ax Go2 xb) 



18 



SGLAZEKRGPEP20120420 printed on January 30, 2013 



+ 



xy 



r. r(123,0) . r(123,0) 
axb J=>t xyb ^* o-yb axy 



X (^01 ax Qoi xy Qo2 yb + Goi ax Qo2 xy Qoi yb + Qo2 ax Goi xy Goi yb) 



axb 

xyz 



■g(124,(234,0)) „(134,(123,0)) 
txyzb ' txyzb 



. r„{124,(234,0)) „(134,(123, 

"T ^tazb \-^taxyz ' -^taxyz 

X Qoi ax Qoi xy Qoi yz 001 zb ■ 



0)) 



+ A 



tayb Bicixy ^tyzb 



(123,0) ^(123,0) 



(B.5) 



The integration yields Eq. (32), which is written in Section 3.3.4 using the 
notation for operators and coefficients that is described below in Appendix 

EH 



Appendix B. 1 Summary of terms and coefficients 
The operators that emerge up to the fourth order in Section [3. 3| are 



Qo2 axb 


= Qoi ax Qoi xb , 


(B.6) 


^03 axb 


= Goi ax Qo2 xb + Qm ax Qoi xb , 


(B.7) 


QqZ axyb 


= Goi ax Goi xy Goi yb i 


(B.8) 


QoA axb 


= Goi ax Go3 xb + G02 ax Qo2 xb + QoZ ax Qoi xb , 


(B.9) 


Soi axyb 


= Goi ax Goi xy Go2 yb + Ooi ax Go2 xy Goi yb 






+ Go2 ax Goi xy Goi yb , 


(B.IO) 


^'04 axyzb 


= QmaxQoixyQoiyzQoizb ■ 


(B.ll) 



The corresponding coefficients are 



„(123,0) _ J^A ^ 

^taxb — / "^^raxfe, 

(124,(234,0)) ^ f dTAraxbB^'^'''\ 
taxyb T axo rxyb ' 

„(134,(123,0)) _ /■* . ,r(123,0) 

^ taxyb ~ Jo '^'^Vb^raxy i 

„(125,(235,(345,0))) _ /"* , 4 „(124,(234,0)) 

axyzb ~ axb -O^ ^y^f. 



„(125,(245,(234,0))) _ / , . „(134,(123,0)) 

axyzb ~ Jo '^"-^^ rxyzb 

„(145,(124,(234,0))) _ /"* . . o(124,(234,0)) 

axyzb ~ Jo "^"-^^ ^axyz 
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„(145, (134,(123,0))) 
t axyzb 



B 



(135,(123,0), (345,0)) 
t axyzb 



[\it A u r(134,(123,0)) 
I u I J-Lj- azb axyz 

Jo 



(It A , r(123,0) o(123,0) 
^rayb ^raxy ^ryzb 



(B.18) 
(B.19) 



The function Afaxb is defined in Eq. (22). The convention for superscripts in 



coefficients B is designed to reflect the origin of corresponding terms. Every 
coefficient has a subscript in the form of a hst of particle configurations. The 
configurations are numbered with natural numbers from left to right. The 
numbers in the superscripts, indicate which configuration appears as a label 
in a corresponding factor under the integral. The first three numbers in a 



superscript refer to subscripts of A defined in Eq. ( 22 ) and the remaining 



numbers refer to the subscripts in coefficients that form other factors under 



the integrals, according to the pattern illustrated by Eqs. (B.12) to (B.19). 
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